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Abstract 

We consider quantum walks on the cycle in the non-stationary case where the 
'coin' operation is allowed to change at each time step. We characterize, in algebraic 
terms, the set of possible state transfers and prove that, as opposed to the stationary 
case, the associate probability distribution may converge to a uniform distribution 
among the nodes of the associated graph. 
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1 Non-stationary quantum walks on the cycle 

Consider a non-oriented graph where all the N nodes have the same degree d and assume 
that, at each time step, a walker makes a choice, out of a set of d elements, {1, ...,d}, a 
(generalized) coin, with probability pi,p 2 , ...,Pd, respectively. The walker starts from a 
given node of the graph and moves in a direction determined by the choice in {l,...,d}. 
After time t, the walker will have a probability P(j,t) of being found in the node j, 
j = 1, N. Such a system is known as a random walk on a graph. A quantum walk is 
the quantum counterpart of a random walk in that both the walker and the coin are seen 
as quantum systems of dimensions N and d, respectively. At each step an operation C is 
performed on the coin system and then an operation is performed on the walker system. 
The latter operation depends on the state of the coin system. 
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Quantum walks have recently received large attention due to the fact that they can 
model quantum algorithms and generate interesting quantum states. There are several 
review papers on quantum walks, their use, dynamics, implementations and generaliza- 
tions (see, e.g., [10], [11]). In most studies presented so far, the coin operation C is fixed 
and repeated at each time step. We shall call this type of quantum walks stationary while 
quantum walks where the coin operation is allowed to change at each time step will be 
named non- stationary. Studies exist on how the parameters of C affect the behavior of 
the quantum walk [19]. The non-stationary case has been considered in both numerical 
and analytic studies where the coin operation is allowed to change at each step according 
to a prescribed sequence or it is random [12], [13], [15], [17]. It is shown that, for cer- 
tain walks, the presence of random noise in C at each step, the so called unitary noise, 
causes a behavior similar to the classical random walk. A non stationary setting can be 
considered also for other types of random walks such as classical walks on groups (as 
for example the Heisenberg group) [7], [8], [14]. These systems are of current interest 
as models of quantum dynamics. The role of the coin process is played by a dynamical 
system which may be characterized by a time varying transformation therefore giving rise 
to a non-stationary random walk on a group. Similar questions to the ones treated here, 
in particular concerning the set of achievable distributions, can be asked in that setting 
as well. 

In this paper, we approach the study of non-stationary quantum walk from the point 
of view of design and control [5]. We consider the coin operation C as a control variable 
which we can change at each step in order to obtain a desired behavior. The first questions 
that arise in this setting are therefore about the type of behavior that can be obtained 
(in particular the type of probability distributions) and whether there are significant 
differences with the stationary case. This paper is a first study in this direction. 

Current proposals for implementations of stationary quantum walks (see [10] and 
references therein) may be modified in order to obtain a non-stationary walk. This is 
discussed for example in [15] for a specific experimental proposal where a variable coin 
operation can be obtained by varying the duration of a laser pulse. 

The quantum walk on the cycle is the simplest finite dimensional quantum walk. The 
study of stationary quantum walks on the cycle was started in [1]. In this case, as a 
consequence of the reversibility of the evolution, the probability distribution P(j,t) does 
not converge to a constant value as t — > oo. This is in contrast with the classical random 
walk on the cycle whose probability distribution converges to a uniform distribution in- 
dependently of the initial state. For this reason, a Cesaro type of alternative probability 
distribution is introduced which is the average of P(j,t) over an interval of time [0, t). 
With this definition, a uniform limit distribution is obtained for number of positions N 
odd, which is independent of the initial state as long as this one is localized in one given 
position. For N even, there is a much richer behavior and different limit distributions are 
obtained for different initial states as discussed in [3], [4]. From an experimental point 
of view, having a uniform Cesaro type limit probability distribution, means that there is 
equal probability of finding the walker in one of the positions by measuring at a random 
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time over a very large interval. 

In dealing with non-stationary walks on the cycle, the first question concerns the 
type of (non-Cesaro) probability distributions that can be obtained. This question is 
of interest in the use of random walks for algorithmic purposes. In fact, there exist 
several computational algorithms which are based on sampling from a given set of objects 
according to a prescribed distribution [18]. These algorithms are referred to as randomized 
algorithms. If one uses a quantum walk to implement one of these algorithms, one can 
obtain the desired sample by measuring the position of the walker. One natural question 
concerns the set of possible distributions available. We shall answer this question in Lie 
algebraic terms in this paper for a non-stationary walk on the cycle and will show that 
it is possible that the probability distribution converges to a uniform distribution. There 
are at least two reasons to consider the uniform distribution with special attention. One 
is that it offers an example of a limit distribution (in the non-Cesaro sense) which is not 
available in the stationary case. In fact we will show that it is possible to reach a separable 
state of the form |1) <g> \ w) where |1) is a state of a two dimensional coin and \w) is a state 
of the walker with all equiprobable positions (cf. formulas (17) and (18) below). Since 
the coin operation is arbitrary, we can set it equal to the identity in the following steps 
and in the following evolution the walker will just move from one position to the other in 
a cycle and the probability of each position will be uniform. The probability distribution 
therefore reaches the uniform value and then stays constant. This is an example of a 
behavior different from the stationary case. The second reason to consider the uniform 
distribution in more detail is that this is the limit of the corresponding classical random 
walk on the cycle. Therefore, if an algorithm uses this feature of the classical counterpart, 
it can be implemented with the non-stationary quantum walk. For example, a randomized 
algorithm which requires sampling from a uniform distribution can be implemented by 
measuring the position of the walker at a large time. If we use the Cesaro definition of 
probability distribution we could perform a measurement but will have to select (again) 
a random time over a large interval. 

A quantum walk on the cycle is a bipartite quantum system C®W, where the system C, 
the coin, is a two level system with ortho normal basis states | + 1) and | — 1). The system 
W, the walker, is an iV-level system with orthonormal basis states |0), |l),...,|iV — 1). 
At the t-th time-step, one performs a coin operation of the form C t <g) 1 where C t is an 
arbitrary (special) unitary operation on the two dimensional Hilbert space associated to 
C, i.e., an element of SU(2). This is followed by a conditional shift S on the Hilbert space 
associated to W defined as 

S\c) ® \j) = |c) <g> \(j + c) mod N). 

By considering the standard basis |e,-), j = l,...,2A r , defined by \ej) := |1) (g> \j — 1), 
and |ej + jv) := | — 1) <8> \j — 1), j = 1, ...,N, the matrix representation of the operator 
C t <S> 1 is C t ®1nxn where Inxn is the N x N identity, 1 ® denotes the Kronecker product 
of matrices, and C t G SU{2). The matrix representation of the operator S is the block 

1 We replace this notation by 1 when there is no ambiguity on the dimension. 
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diagonal matrix diag (F,F T ), where F is the basic circulant permutation matrix, that 
is, 

1\ 




/0 
1 
1 



S:= diag(F,F T ), 







(1) 



Vo • • • • i o/ 

The probability of finding the walker in state \j — 1), j = 1, ...,N, is the sum of the 
probabilities of finding the state of the composite system W in |1) <g) \j — 1) := |1, j — 1) 
and | — 1) <g) \j — 1) := | — 1, j — 1). That is, if is the state of the composite system, 

P(j -l,t) = - 1>| 2 + \m)\ - U - l)! 2 = I W)h>| 2 + I W*)te + *>| 2 . (2) 

By writing 

IV>) : = S a fcl e fc)' H l^fcl 2 = !> ( 3 ) 

fc=l fc=l 

we have 

P(j - 1, t) = \ aj \ 2 + \a 0+N \\ j — 1, ...,N. (4) 
In what follows we shall make the following standing assumption. 

Assumption: N is an odd number. 

We use this technical assumption in several steps and in particular in Theorem 1 to 
show that the matrix S defined in (1) is in a certain Lie group (cf. formula (11)). In the 
stationary case, problems with the iV even case arise from the degeneracy of eigenvalues 
in the basic S(C ® 1) operation. As we have mentioned, this leads to a very different and 
richer behavior with respect the iV odd case. 



2 Characterization of the admissible evolutions 

In this section, we characterize the set of unitary evolutions available for a non-stationary 
quantum walk on the cycle, that is, the set of available state transfers. This is the set 
of finite products of operators of the form S(C <S> ljvxiv) where S is defined in (1) and 
C G SU(2). We denote such a set by G. The set G is a group. It is in fact a Lie group 
as it is shown in the following theorem. In order to state this theorem, we need to recall 
some properties of circulant matrices [6] and define two Lie algebras. 

Circulant N x N matrices with complex entries form a vector space over the real 
numbers. Each matrix is determined by the first row since all the other rows can be 
obtained by cyclic permutation of the first one. Moreover every complex circulant matrix 
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R can be written as linear combination with complex coefficients of the basic permutation 
matrix F defined in (1) and its powers from to N — 1, i.e., 



JV-l 

i?:=E«^> (5) 

1=0 

with N complex coefficients a$, Qjv-i- All the circulant matrices commute. If we require 
that R is not only circulant but also skew- Hermit ian then we must have 

JV-l N-l 

Rl = a*l + E < FlT = ~R = -«ol - E a i pl i ( 6 ) 
i=i i=i 

and with a change of index / — > N — I and using F N ~ lT — F l , we have 

JV-l JV-l 

= a*l + E a *N-iF l = -oqI - E aiF 1 - (7) 
i=i i=i 

This gives the relations 

N-l 

a* = -a , a * N _ l = -a h 1 = 1,..., — - — . (8) 

Equations (7) constitute iV independent relations on the 2N real parameters of R and show 
that the space of skew-Hermitian circulant matrices is a real vector space of dimension 
N. 

Denote by C the Lie algebra spanned by the 2N x 2N skew-Hermitian matrices of the 
form 

R o \ _ A T / Q 



Li:= U -Rj and L2:= l-Qt oj' (9) 

with i? a skew-Hermitian circulant N x N matrix and Q a general circulant matrix. It 
is easily seen that this is in fact a Lie algebra of (real) dimension 3iV; the fact that it is 
closed under Lie bracket being a consequence of the fact that the product of two circulant 
matrices is another circulant matrix. Notice, in particular, that matrices of the type L\ 
form an Abelian subalgebra of dimension N. We denote by e c the connected Lie group 
associated to C. 



Theorem 1 The set G of possible evolutions of a non- stationary quantum walk is the 
Lie group e c . 

Proof. We define an auxiliary Lie algebra CJ , prove that G = e c ' and then prove that 
C = £'. The claim then follows from the correspondence between Lie groups and Lie 
algebras. We denote by CJ the Lie algebra generated by the set 

T := {su(2) ® 1, S (su(2) 1) S T , S N ' X (su(2) ® 1) S (N ' 1)T }, (10) 
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where S is defined in (1) and T denotes transposition. 

To show that G C e £ ', it is enough to show that both C <g> 1, C G SU (2), and S are in 
e c . This fact is obvious for C<g> 1, since this is the exponential of an element in su{2) <g> 1. 

For S, we consider the elements ^ ^ ^ <E> 1 and S 1 ^ - - (( , „ ] (g) 1 ] S^) 7 , both 

in e £ ', and calculate with (1) 



-1 

1 



(ii) 



, F (N-1)T \_{F 

I F^" 1 J ~ I F T 



:= 5. 



We have used F^ N ~^ T = F. 

To show that e c ' C G it is enough to show that every element of the type S j (X <g> 1) S^ 7 , 
with X G SU(2), j = 0, ...,iV — 1, can be written as the finite product of elements of 
the form S (C <g> 1) with C G S'[/(2). 2 This is readily seen because, with X G 577(2), for 
every j, 

S j (X l NxN ) S lT = (S j (X l NxN )) (s N ' j (l 2x2 ® ljvxiv)) . (12) 

To conclude the proof, we show that C = £ showing that JF C £ and a basis of C 
can be obtained as (repeated) Lie brackets and-or linear combinations of elements of T 
in (10). A general matrix in T has the form, with A G su(2), 

S-* (A <S> Inxn) S^ t = (13) 

F j \ / ibl NxN al NxN \(Fi T 0\ = f ibl NxN aF 2 i 
F jT ) \—a*l NxN -ibl NxN )\0 F^J \ —a*(F 2jT ) -ibl NxN 

with general b real and a complex, j = 0, N — 1. This is clearly in £. Elements of 

jF k 



the form L 2 in (9) are real linear combinations of elements of the form ^ ^*pkT 
which are of the form in (13) with b = 0, 7 = a and j : = | for A; even and j = for A; 



odd. A basis for the real elements of the type L x is given by the linearly independent 
elements 



F j ~F jT \ ._ N-l 

—(F j — F jT ) J ' J" 1 '-' 2 ' 1 j 

IatxaA , / F J 



These are obtained as Lie brackets of ) and T which are 

\— IatxAT / \—r J (J / 

both of type L 2 . A basis for the purely imaginary elements of type L\ is given by the 
linearly independent elements of the type 



(F J + F jT ) \ . _ N 



+ F^' T ) 



j = 0, , (15) 



2 Recall that every element of a connected Lie group can be obtained as the finite product of expo- 
nentials of a set of generators of the corresponding Lie algebra (see, e.g., [9]) and the exponential map is 
surjective on SU(2) (see, e.g., [16]). 
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which are obtained as Lie brackets of 
pletes the proof of the theorem. 











and 





Hnxn 



HnxN 





This com- 
□ 



3 Obtaining the uniform distribution 

The Lie group G = e c , having dimension 37V, is not isomorphic, for N > 3, to SU(2N) 
(which has dimension AN 2 — 1) nor to Sp(N) (which has dimension N(2N + 1)). Therefore 
G = e c is not transitive on the complex sphere of dimension 2N which means that there 
are state transfers for the quantum system of coin and walker which are not induced by 
any transformation in G [2]. Some state transfers are of special interest. In particular, 
we are interested in whether a state of the form 



\lpin) ■= \lpcoin) ® |0), 



(16) 



that is, a state corresponding to the walker with certainty in position |0), can be trans- 
ferred to a state corresponding to the uniform distribution. This is a state where the 
probability P(j — l,t) in (2) is equal to j^, for every j = 1,...,JV, at some t, that is, 
the walker is found in any position with the same probability. Since, VC G SU(2), 
C <S> Inxn ^ we can assume, without loss of generality, that \ip CO i n ) m (16) is |1) so 
that the problem is to transfer the state |ei) := [1,0, ...,0] T to a state with the desired 
property. We shall show in the following that such a state transfer is possible. 

Theorem 2 There exists a matrix L in C such that 

(n\ 

T2 



lei) 



r N 



V o J 



;i7) 



where 



l r i| 



\ r 2\ 



\ r N\ = T7 



1 

N' 



(18) 



In order to prove this theorem we first prove a lemma. Recall the definition of the 
Fourier matrix $ of order iV (see, e.g., [6]). This is defined so that its conjugate transposed 
is 





fl 


1 


1 


1 


1 


\ 




1 




uj 2 


a; 3 


u N - 1 




1 


1 


CO 2 


cu 4 


u 6 








1 


u 3 






W 3(JV-1) 






u 






W 3(JV-1) _ 


JN-1KN- 





(19) 
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where u is the N-th root of the unity, that is u := e % n . The Fourier matrix $ is unitary. 
Lemma 3.1 Define 



xi : = 



1(1 - 1)tt 
~~N ' 



l = 0,l,...,N-l. 



Then 



/ri\ 




/ e « \ 


r 2 








:= > 






\/iV 




\r N ) 







has the property (18). 

Proof. From (19) and (21), we obtain 



N-l 



This, using the definition of u>, gives 



UJ 



(h-l)l ix t 



h = l,...,N. 



N-l 



1 / - i2*(h-l)l 

Th = n \ n e " e 



We calculate |r/j| 2 , h — 1, ...,iV, as 



-, AT-l 



AT2 



E - 

h,J2=0 



AT-l 



(20) 



(21) 



(22) 



(23) 



(24) 



J_ _|_ A ^ R e (e^Ca-JiXft-lJgttea-x,^ _ 

^ ^ 2 {h,h}=o 

The sum in the last term is intended over all the pairs of indices {h,l 2 }, with li ^ l 2 , 
where only one is chosen between {7 1; l 2 } and {l 2 , h}. Because of the presence of the real 
part ' Re ' it is not important which pair is chosen. We now show that, with the choice 
(20), the last term of this expression is zero for every h, which will prove the claim that 

M 2 = jr- 

It is convenient to re-write the sum by regrouping elements corresponding to l 2 — l\ = 
p mod N, for p = 1, N — 1. This means l 2 — l\ = p or l\ — l 2 = N — p. We have 



AT-l 



53 Re (e i ^ (,a ~ ,l)(h ~ 1) e i(x, 2~ x, i ) ) 

{h,h}=0 



(25) 



AT-l 



V Re V e i(i2-h)(h-i)f e i(x h -x h ) + e i(i 2 - h ) {h -i)f e i( Xl2 - Xh ) 



P =i 



\h—h=p 



h-l 2 =N—p 
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Doing the substitution li = I and I2 = l+p in the first term of the sum and the substitution 
li = I and l 2 = I — (N — p) in the second term, this sum becomes 



N-l 



'N-l-p 



N-l 



(26) 



P =i 



1=0 



=N-p 



We now show that, with the choice (20), the content of the innermost parenthesis in the 
above expression, i.e., 



N-l-p 



N-l 



M := M(p) := Y e i( - Xl+p ~ Xl) + Yl e^ Xl -<- N -^~ Xl \ 



(27) 



1=0 



l=N- 



is zero for each p which will conclude the proof of the Lemma. Replacing (20) in (27) and 
after some algebraic manipulations, we obtain 



»t/ \ N x^ P i'^( p(p -^ 1 pi) , ji* ( ( N -p)<- N -p+v 

M(p) = Y 6 ^ 2 ' + 2^ e 2 

1=0 l=N-p 



N-l 

M^e^(^). (28) 
1=0 



Thus we have 



2tv p(p-l) x ■ 2irpl -2__ 

I N 2 



1=0 



7t p(p-i) 1 -e i2np 



1 — e n 



Wp ^ mod AT. (29) 



This concludes the proof of the Lemma. 



□ 



We are now ready to prove Theorem 2 
Proof. (Proof of Theorem 2) We choose L as a matrix of the form L x in (9) so that e L 
has the form 

(30) 



e R 
e 



-R 1 



with a general skew-Hermitian N x N circulant matrix. The problem is therefore to 
find a circulant matrix R so that 
















e R 










W 



(31) 



with ri,...,rN satisfying (18). Any circulant matrix R is diagonalized by the Fourier 
matrix (19) of the corresponding dimension, that is, 



R = $ f A$, 



(32) 



with A diagonal. Conversely every matrix of the form on the right hand side is circulant 
[6]. Moreover if A = diag (i\ , zAjv-i), with A;, / = 0, AT — 1 real numbers, R is 
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skew- Hermit ian. In this case, we have 



R 























= &e A l_ 


i 

















Ve 



(33) 



/ 



Choosing A; = xi, I = 0, ...,N — 1 with the definition (20), the theorem follows from 
Lemma 3.1. □ 

Other states with the same property can be obtained by applying a transformation 
U <g> 1, U e SU(2), which is in G. In particular notice that the state (17) is a separable 
state. 



4 Conclusion 

Non-stationary quantum walks have properties which distinguish them from stationary 
ones. Moreover they are amenable of study with the methods of quantum control. In fact, 
several problems, such as obtaining a given evolution, can be seen as control problems 
where the evolution of the coin plays the role of the control. In this paper we have shown 
that, opposite to the stationary case, a non- stationary quantum walk on the cycle may 
converge to a constant distribution and in particular to a uniform distribution as for 
classical random walks. A constructive approach to achieve this and other evolutions of 
interest for general quantum walks will be the subject of future research. 

Acknowledgement D.D. is grateful to Mark Hillery for helpful conversations on the 
topic of quantum walks. The authors also thank the reviewers for useful suggestions. D. 
D'Alessandro's research was supported by NSF under Career Grant ECS-0237925. 
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